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We show the existence of a rather general class of closed cosmological models of Bianchi type IX 
that do not exhibit recollapse but expand for all times. This is despite the fact that these models 
satisfy the strong energy condition by a wide margin. 
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In relativistic cosmology, the trichotomy of the Fried- 
mann-Robertson- Walker (FRW) models is of prime im- 
portance. The spatial geometry determines the evolution 
of the cosmological model: In the hyperboloidal (fc = —1) 
or the spatially fiat (fc = 0) case, the universe exhibits 
an initial singularity ('big bang'), and from that 'mo- 
ment' on, the universe is forever expanding. In the case 
of a closed cosmological model (fc = +1) we observe a 
fundamentally different behavior: "[...] the dynamical 
equations of general relativity show that the spatially 
closed 3-sphere universe will exist for only a finite span 
of time. [. . . ] at a finite time after the big bang, the 
universe will achieve a maximum size [...], and then will 
begin to recontract. [. . . ] a finite time after recontrac- 
tion begins, a 'big crunch' will occur." The quotation is 
taken from 12 . The recollapse of closed FRW cosmolo- 



gies holds because the strong energy condition is imposed 
on the matter (which is assumed to be a perfect fluid), 
i.e., p + 3p > 0, where p is the energy density and p the 
pressure. 

On this basis it would be tempting to view the recol- 
lapse of the spatially closed FRW cosmologies satisfying 
the strong energy condition as a paradigm for models 
with less symmetries. That this belief is erroneous has 
been demonstrated in At least in the case of cosmo- 
logical models of Bianchi type IX, into which the closed 
FRW models are naturally embedded, there exists a rig- 
orous result by Lin and Wald 8]: Assuming the domi- 
nant energy condition, i.e., \pi\ < p for the anisotropic 
pressures p±, P2, pa, and non- negative average principal 
pressure p, i.e., 3p = (pi +P2 + £3) > 0, then the 'closed- 
universe-recollapse conjecture' [3| holds. In the locally 
rotationally symmetric (LRS) case with isotropic matter, 
it is sufficient to require that p + 3p > ep for an arbi- 
trarily small e > 0, see Q- (Note that if p/p — » —1/3, 
there exist models that expand forever approaching the 
Einstein static universe in the limit.) 

In this letter, we approach the problem from a different 
direction by investigating cases where 'closed-universe- 
recollapse' does not hold. We prove that there exist 
cosmological models (of Bianchi type IX) satisfying the 
strong energy condition that do not recollapse but ex- 
pand forever. Two points are important to emphasize: 



(i) For these models the strong energy condition is sat- 
isfied 'by a wide margin'. The assumption we make is 
that uu = p/p is a constant, i.e., w — const > —1/3; 
however, w < (1 — y/3)/3 s» —0.244, hence the average 
pressure is not positive, (ii) We prove the existence of a 
typical class of models that expand forever (where 'typ- 
ical' is understood in the sense of an open set of initial 
data of the Einstein equations). An interesting obser- 
vation is that these cosmological models exhibit partial 
(i.e., directional) accelerated expansion for late times. 

As a matter of course, we do not propose the cosmolog- 
ical models we analyze as actual models of the universe. 
However, we want to em pha size that the matter model 
we consider is not 'exotic' [10] . i.e., it satisfies all the stan- 
dard energy conditions (weak, strong and dominant), as 
opposed to certain 'exotic' matter models in cosmology 
(e.g., phantom fields and dark energy, see for instance 
[l|, |9[ ; the breaking of the energy conditions stems from 
the aim to account for the accelerated expansion of the 
universe). The properties of the matter source we con- 
sider in this paper resemble those of collisionless mat- 
ter, elastic matter, and magnetic fields as regards their 
fundamental aspects [1, Hj]. Classes of models encom- 
passing these important examples (or a subset thereof) 
have been the basis of previous work, see, e.g., [ij. An- 
other explicit example, an example that satisfies our con- 
crete assumptions, is the anisotropic fluid model Q. For 
the models we consider the anisotropic pressures (par- 
allel and perpendicular pressure) are required to satisfy 
certain bounds (that are compatible with the energy con- 
ditions). In particular, the isotropic pressure (average 
pressure) depends linearly on the energy density, as is 
usual for perfect fluids, where the proportionality con- 
stant is strictly larger than — |. Note, however, that the 
approach we take does not require to specify a concrete 
matter model as long as the basic assumptions of 0, Q 
and the necessary bounds on the anisotropic pressures 
are satisfied. Finally, note that the restriction to matter 
sources of the general type of [!, [j| might yield a special 
case of the 'closed-universe-recollapse conjecture' and, in 
view of the results of this paper, lead to specific bounds 
on the matter quantities. 

Let us briefly comment on the approach we take. We 
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use the dynamical systems approach to spatially ho- 
mogeneous cosmologies 11 1. However, as we will see, 



it is essential to avoid the standard Hubblc-normalized 
variables — the results we present here are rather elusive 
in that approach. 

Equations. For models of Bianchi type IX the metric 
can be written as 



ds 2 = -dt 2 + gij(t) 



(1) 



where {a) 1 , cj 2 , d) 3 } is a symmetry-adapted coframe sat- 
isfying dili 1 = —uj 2 A Co 3 (and cyclic permutations). The 
Einstein equations comprise evolution equations for the 
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and for the extrinsic curva- 
The Gauss constraint reads 



metric, <9t<7ij 
ture fc 4 -, see, e.g. 
3 R + (tr k) 2 — k % ^k? i = 2p , where p = — T° is the energy 
density associated with the energy-momentum tensor T^ v 
and 3 i? the scalar three-curvature. In the vacuum case 
or orthogonal perfect fluid case the metric is diagonal, 
i.e., gij(t) = diag (ffii(t),322(*),ff33(*)); in the locally ro- 
tationally symmetric (LRS) case we have g>22 (t) = 93i(t)- 
We restrict ourselves to diagonal metrics even if the mat- 
ter source is anisotropic. 

In the diagonal case, an isotropic matter source is char- 
acterized by an energy-momentum tensor with T l j — 
diag \Pi,P2,P3)- The 'isotropic pressure' p is the aver- 
age of the anisotropic pressures pi, p%, ps, i.e., trT = 3p. 
Define w and Wi, i = 1, 2, 3, according to 



p = wp. 



Pi 



Wip ; 



(2) 



obviously, w\ + u>2 + VJ3 = 3u>. Matter that is consistent 
with LRS symmetry satisfies W2 = 103. For perfect fluids, 
w\ = 11)2 — ^3 = w, where w is typically assumed to be 
a constant. In this paper we consider anisotropic matter 
that generalizes perfect fluid matter: We assume that 
the energy density and the isotropic pressure satisfy a 
linear equation of state, w — const, where the strong 
energy condition is supposed to hold, i.e., w > —1/3. The 
rescaled anisotropic pressures, wi, u>2, W3, are assumed 
to be functions of the metric via (s±, S2, S3), where 



Sfc 



g kk (g 11 + g 22 + g 33 ) 1 (no sum over k) ; (3) 



obviously, Si + S2 + S3 = 1. The functions 

W k = Wk(s 1 ,B2,8a) 



(4) 



are such that there exists an isotropic state of the mat- 
ter where w\ = W2 = w$ = w, and remain bounded 
(and take limits) under extreme conditions (when one 
or more of the Sj are zero). In particular, there exists 
a constant t>_ such that wi(0, S2, S3) = 102(^1,0, S3) = 
1^3(31,82,0) — V-. There exist excellent examples for 
matter models of this type, e.g., collisionless matter, elas- 
tic matter, or magnetic fields; for a detailed discussion we 
refer to lal. 



Let us define the Hubble scalar H as H = — trfc/3 
and the shear tensor (ai, 02, 03) as the traceless part of 
the extrinsic curvature, i.e., k l i = —H — ai (no sum over 
i); o\ + (72 + 03 = 0. Furthermore, we introduce the 
'densitized metric' (n\, ri2, n.3) by rik — gfcfc(det <?) -1 ^ 2 ; 
note that > 0. Then the Einstein equations can be 
expressed as evolution equations for H, (01,02,03), and 
(m, ^21^3) plus one constraint, which can be used to 
express p in terms of the other variables; this leads to 
the fact that the matter enters the equations only via w 
and Wi, W2, W3. 

In the LRS case, which we will focus on henceforth, 
there exists a plane of rotational symmetry, which we 
choose to be spanned by the second and the third frame 
vector. Accordingly, ri2 — TI3 and 02 = 03 as well as 
S2 = S3; consistently, the matter satisfies W2 = W3. Let 
o := 02 and let s :— S2; then 

(01, 02,0-3) = (—2a, a, a), (si, s 2 , s 3 ) = (l-2s, s, s). (5) 

Eq. © implies s = (2 + n 2 /n 1 )- 1 , so that s G (0, 1/2). 
Finally we abbreviate the rescaled anisotropic pressure in 
the plane of rotational symmetry by u; more specifically, 



u(s) := 102(1 ~ 2s, s, s) 



(6) 



The Einstein equations can then be expressed in the vari- 
ables H, o, ni, and ri2, where u(s) appears in these equa- 
tions. 

In the dynamical system approach to cosmology the 
Einstein equations are expressed in terms of normalized 
variables. We define the 'dominant variable' D, see, 
e.g., 0, by 



D = JH 



nin 2 



(7a) 



and we introduce normalized variables according to 
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in addition we use the variable 



( 2 +;r 



(7b) 



(7c) 



Further, we define a normalized energy density Qjj = 
p/(3D 2 ), and we replace the cosmological time t by a 
rescaled time variable r through 



d _ 1 d 
dr ~ ~D~dt' 

Henceforth, a prime denotes differentiation w.r.t. r. 
Remark. The evolution equation for H is 

H' = -^(H 2 + q D D 2 ), 



(8) 



(9) 
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where q D is given by q D = 2E 2 , + (1/2) (1 + 3w)fl D . This 
leads to an important remark: Equation ([9]) implies that 
H is decreasing if H — 0, i.e., H'\ H=0 = —q D D~ x < 0; 
therefore, a cosmological model with H(t ) > at some 
time t satisfies H(t) > Vr < tq. Consequently, by 
proving the existence of models that satisfy H (r) > for 
all sufficiently large t, we prove the existence of models 
with H > and thus positive expansion for all times. 

It is not difficult to prove that the transformation ([7]) 
between the 'metric variables' H, a, m, n 2 and the dy- 
namical systems variables D, E^, r, s is one-to-one on 
the set H > 0. (This is sufficient for our purposes, see 
the previous remark.) 

Expressed in the variables D, So, r, s, the Einstein 
evolution equations split into a decoupled equation for 
D, 



D' = -D(H D (1 + q D ) + Ed(1 - H 



D) i 



(10) 



and a system of coupled equations for the normalized 
variables ([7b| and (l7cl). 



54E D s 2 



(11a) 



r' = r{2H DiqD -H^ D )- % ^ + ^ ; , 

E' D = -(2 - q D )H D Z D - (1 - - E|,) + 

+ §JV?d + 3^ («(«)-«>), (lib) 



s' = -12s (h ~ alS ; 



(11c) 



where g D = 2EJ, + ±(1 + 3w)fi D , and i? D = H D {r,s) 
and _/Vi£> = N\d{t, s) are functions of r and s, see (|T3"|). 
The Gauss constraint reads 



K + MiD(r,s) + n D = l 



(12) 



it is used to solve for fi^. The functions H]j(7\s) and 
N 1D (r,s) are 



H D :=^- = H D (r,s) = Jl-r 1 2 " 



VI - 4s + 13s 2 



(13a) 



Nw:=ji=N 1D (r,s) = J 71 fe= 2 ; (13b) 

y Vl - 4s + 13s 2 

in particular, Hd(t, s) and Nuj(r, s) are well-defined and 
regular on the preimage of the set R + x R + . (We refrain 
from going into details in this paper, since we merely use 
that (fl~3"|) is well-behaved for sufficiently small r; however, 
we may refer to Q.) 

Remark. More common than ([7]) is the Hubble-norma- 
lized approach, see, e.g.,[ll|, where the Hubble scalar 
H is employed instead of D to construct scale-invariant 
variables and a simpler set of normalized variables is used 



instead of (l7b|) and ((7c)) . Although the resulting equa- 
tions are simpler than (|11[) . we do not have a choice; we 
will see that the Hubble-normalized approach necessarily 
fails to uncover our results. 

The dynamical system (fTTj) completely describes the 
dynamics of locally rotationally symmetric cosmologi- 
cal models of Bianchi type IX in their expanding phase 
(H > 0). In other words, each solution of (fTTj) yields an 
LRS Bianchi type IX model in its expanding phase, and 
conversely, the expanding phase of each model is repre- 
sented by a solution of (fTTj) . The main advantage of the 
system (fTT]) over other representations of the Einstein 
equations lies in its extendibitly: The system (fTTj) pos- 
sesses a regular extension to its boundaries Ed = ±1, 
s = 0, s = i, and r = 0. 



Anisotropic matter. The function u(s) in (|llb[) en- 
codes the properties of the matter model. For isotropic 
matter we have u(s) = w; for anisotropic matter, u(s) 
represents the rescaled anisotropic pressure in the plane 
of local rotational symmetry, cf. ©; recall that s £ 
[0,1/2]. From u>2(si, 0, S3) = w_ we obtain u(0) = 
v-, cf. © and ©. The value of u at s = 1/2 
is not independent; to see this recall first that w\ + 
2w2 = iw, now, s = 1/2 corresponds to (si,S2,S3) = 
(0, 1/2, 1/2), so 101(0, 1/2, 1/2) + 2w 2 (0, 1/2, 1/2) = 3w; 
hence, wi(0, S3) = V- results in t>_ + 2u(l/2) = 3w. 
Summarizing, 



u(0) 



u(l/2) 



w 



(1/2) (t 



(14) 



For reasonable matter models like collisionless matter, 
elastic matter, or magnetic fields, u(s) is a monotone 
function on [0, 1/2] interpolating between the values (fT4| : 
we refer to H for examples. It is often beneficial to use 
the anisotropy parameter (3 instead of the constant v_; 
it is defined as 



2(w - v-) 



1 



w 



(15) 



Finally, let us state the assumptions on the matter 
model we consider in this paper. We assume that 

~l < w < ^—^^ -0.244, (16) 

hence the strong energy condition is satisfied. Further- 
more, V— is assumed to satisfy 



i (l + 6w - V-3 + (1 - 3w) 2 ) < v_ 

v- < i (l + 6w + V-3+ (1 - 3w) 2 ) . 



(17) 



The admissible domain of the parameters w and w_ is 
depicted in Fig. [TJ Clearly, the dominant energy condi- 
tion is satisfied since the rescaled anisotropic pressure V- 



4 



v± 





1 

3 


-0.30 


-0.27 3 


W 










- -0.08 








**> 


l 

- -0.58 




v + 






- -1 



FIG. 1: The admissible values of the rescaled isotropic pres- 
sure w and V-, v+ = 3w — 2v- , which represent (the extremes 
of) the rescaled anisotropic pressures in the plane of symme- 
try and orthogonal to it, respectively. Both the strong and the 
dominant energy condition are satisfied. (1 — v3)/3 ~ —0.24. 



(which is the pressure in the plane of symmetry) and its 
counterpart v + = 3w — 2v- (which is the pressure in the 
orthogonal direction) satisfy |v±| < 1. 

Results. Since the dynamical system (|TT|) extends 
regularly to r = it is suggestive to analyze the system 
induced on that surface. From (|13[) we obtain Ho\ r =o = 
1, Nf D \ r=0 = 0, so that O becomes Q D = 1-Y? D . This 
in turn implies qo = |(1 + 3u>) + |(1 — w)Y? D . Insertion 
into (HU) yields 

T,' D = -3(1 - Y? D ) - w)^ D - (u{s) -vifj, (18a) 
s' = -12Srjs(i-s) . (18b) 

The state space for this two-dimensional dynamical sys- 
tem is [—1,1] x [0, 1/2]. The dynamical systems analysis 
is straightforward. We use that u(s) is a function such 
that u(0) = V- and u(l/2) = w + (w — v-)/2; here, w and 
V- are assumed to satisfy (fT6|) and (fT7|) . respectively. 

We focus our attention on the fixed point R of ([18)1. 
which is given by 

R : r = 0, s = 0, Sjj = -ft . (19) 

It is straightforward to prove that R is a sink for the flow 
of the system (fl8|) , because 

s"V| R = 6/3 < 0, (20a) 
(E Z5 +/3)- 1 (E D +/3)'| R = -|(l-/3 2 )(l- w ) < 0, 

hence the eigenvalues of the linearization of (fT8|) at R are 
negative. 

The crucial property of the fixed point R is revealed 
by considering the full system ([lip : R is a sink not only 



on the boundary r = 0, but also for the full system (jlll) . 
To see this we simply compute 

r~ V| R = 2(§ (1 - w)/3 2 + |(1 + 3w) + < (20c) 

and use (fl7|) to establish that the r.h. side is negative. 

Since the fixed point R is a sink for the flow of the 
system ([TT[) . i.e., for LRS Bianchi type IX models with 
anisotropic matter satisfying (fT6|) and (TT7| . there exists 
an open subset of LRS type IX initial data such that the 
corresponding solutions converge to R. Because H D = 1 
at R we infer from ([13ap that H is positive for these 
solutions for late times; by the remark following © we 
obtain that H is positive, i.e., these models are expanding 
for all times. 

In the following we analyze the asymptotic behavior of 
these forever expanding LRS type IX solutions in terms 
of the metric variables: For a solution of ([111) that con- 
verges to the point R we find r — > 0, s — > 0, Su — > —/J, 
hence Hd — > 1, Nm — > by IT3|) . Furthermore, 
N2D '■= ni/D — > 00; to see this we first note that s — > 
implies s ~ NidN^jj by (f7c[) . Then, from iViu ~ y 7 ^, 
see (I13bp . we conclude that N 2 d ~ y/r/s. Using that 
s~V|r = 6/3, cf. ([2TE[I . and r-V| R = 2(g + /3), cf. ([20c) . 
where 

?:=g2 3 | R =|(l-«;) j 8 2 + |(H-3u;), 

we finally get ^V^^dIr = 9^2/3 > 0, and the claim 
follows. The fact that N2D = n^/D — > 00 implies that 
N2 := TI2/H — > 00 as r — > 00 since H c± D m the 
limit (because = iJ/Z) = 1). This property is com- 
pletely unproblematic here, but makes the treatment of 
the problem extremely difficult in the standard Hubblc- 
normalized approach. 

For large r, the shear variables are a\ = —2a = 
-2DY, D ~ 2D/3 and a 2 = a z = a = DY> D ~ -L>/3. 
To obtain the metric we integrate dtgkk = ^gkk{H + 0"^) 
(which is a consequence of c>t(?ij = —2guk l j). In the first 
step we note that H = HrjD = D so that the equation 
reads d t gkk = ^9kk{D + a k ), i.e., 

6^11 = 2311(1 + 2,8), ^22 = 2322(1-/3). (21) 

Second, we integrate _D _1 Z)'|r, = — (1 + q) and ([5]) to get 
t = (1 + g) -1 log(< — tg) + T a with constants ro and tg. 
Therefore, (12T|) leads to an asymptotic behavior of the 
metric represented by 

2(l + 2ff) 2(1-3) 

3ii oc (t - t ) i+« , 322 oc (t - t ) !+« (22) 

as t — > 00. Note that /3 is negative, but 1 + 2/3 is posi- 
tive, which is a consequence of (jT7j) . An interesting ob- 
servation is the occurence of partial (directional) accel- 
erated expansion. A straightforward calculation reveals 
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that (1— ft)/(l+q) > 1, hence lengths in the plane of local 
rotational symmetry expand at an accelerated rate. The 
maximal rate of acceleration is obtained by maximizing 
(1 — /3)/(l + q) over the domain depicted in Fig. [TJ the 
maximal value of (1 — /?)/(l + q) « 1.112 is attained for 
w close to —1/3. 

We conclude by restating the main result: The behav- 
ior (|22|) is typical, i.e., there exists an open set of LRS 
type IX initial data such that the associated solutions 
of the Einstein equations with anisotropic matter behave 
as (|22p as t —> oo and thus expand forever. (These so- 
lutions do not behave extraordinarily in other respects; 
for instance, there exist forever expanding solutions that 
isotropize toward the singularity.) In this sense, eternal 
expansion is as likely as recollapse in the case of LRS 
Bianchi type IX with anisotropic matter that satisfies 
the conditions ([TBI) an d (fl7|l- and thus in particular the 
strong energy condition. 
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